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ABSTRACT 

We present an initial sample of 19 intermediate-mass black hole candidates in active galactic nuclei culled 
from the First Data Release of the Sloan Digital Sky Survey. Using the linewidth-luminosity-mass scaling 
relation established for broad-line active nuclei, we estimate black hole masses in the range of Mbh ~ 8 x 
10 4 - 10 6 M Q , a regime in which only two objects are currently known. The absolute magnitudes are faint for 
active galactic nuclei, ranging fromM g «— 15 to —18 mag, while the bolometric luminosities are all close to the 
Eddington limit. The entire sample formally satisfies the linewidth criterion for so-called narrow-line Seyfert 1 
galaxies; however, they display a wider range of Fe II and [O III] A5007 line strengths than is typically observed 
in this class of objects. Although the available imaging data are of insufficient quality to ascertain the detailed 
morphologies of the host galaxies, it is likely that the majority of the hosts are relatively late-type systems. The 
host galaxies have estimated g-band luminosities ~ 1 mag fainter than M* for the general galaxy population 
at z ~ 0.1. Beyond simply extending the known mass range of central black holes in galactic nuclei, these 
objects provide unique observational constraints on the progenitors of supermassive black holes. They are also 
expected to contribute significantly to the integrated signal for future gravitational wave experiments. 
Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert 



1. INTRODUCTION 

Dynamical studies have established the existence of super- 
massive black holes (BHs) with masses M BH « 1O 6 -1O 9 M 
in the centers of most, if not all, local galaxies with bulges 
(Magorrian et al. 1998; Richstone 2004). A significant chal- 
lenge to any model of cosmological BH growth is the na- 
ture of seed BHs. Observations of intermediate-mass BHs 
(Mbh ~ 1 3 - 1 6 M Q ) in the local Universe would provide the 
most direct empirical constraints on such seeds. The Mbh-o* 
relation (Ferrarese & Merritt 2000; Gebhardt et al. 2000a; 
Tremaine et al. 2002) strongly indicates feedback in the co- 
evolution of galaxies and BHs, but it is still an open ques- 
tion whether the relation is established at early or late times. 
If, as some models suggest (e.g., Di Matteo et al. 2003), 
the Mbh _ c+ relation is established late in the evolution of 
a galaxy, then we may expect to see deviations from this re- 
lation in a population of intermediate-mass BHs that is not 
yet fully grown. There are other motivations to characterize 
the lowest end of the nuclear BH mass distribution. Certainly 
their number and mass distribution will impact the expected 
integrated background in gravitational radiation (e.g., Hughes 
2002). Presumably they will pass through an active phase of 
accretion and contribute at some level to the observed inte- 
grated background radiation, particularly at X-ray energies. 
The range of accretion phenomena they display will depend 
on the major accretion modes for BHs of such low mass. The 
accretion states and, by extension, radiative growth mecha- 
nisms of this population will therefore have direct implica- 
tions for the growth history of primordial BHs. Apart from 
their physical implications, intermediate-mass BHs have prac- 
tical value in that they offer tremendous leverage in anchoring 
local BH-galaxy scaling correlations such as the Mbh _ re- 
lation. 

Recently discovered intermediate-mass BH candidates in 



galactic nuclei offer tantalizing hints as to the nature of this 
population. NGC 4395 is a very late-type (Sdm) spiral with 
no bulge component, whose central stellar velocity dispersion 
is ov< 30 km s -1 (Filippenko & Ho 2003). Nevertheless, it 
has the emission properties of a type 1 active galactic nu- 
cleus (AGN), with broad permitted optical and UV emission 
lines (Filippenko & Sargent 1989; Filippenko, Ho, & Sargent 
1993), a compact radio core (Ho & Ulvestad 2001) with a 
nonthermal brightness temperature (Wrobel, Fassnacht, & Ho 
2001), and a pointlike, hard X-ray source (Ho et al. 2001) that 
is highly variable (Moran et al. 1999, 2004; Shih, Iwasawa, 
& Fabian 2003). Mass estimates based on the H/3 linewidth- 
luminosity-mass scaling relation and X-ray variability suggest 
a BH mass of 10 4 - 10 5 M (Filippenko & Ho 2003). In- 
terestingly, these agree with the limit of < 10 5 M Q derived 
from the M BH - cr* relation (Filippenko & Ho 2003). POX 
52, on the other hand, has a dwarf elliptical host. Recently 
revisited by Barth et al. (2004), POX 52 was first observed 
by Kunth, Sargent, & Bothun (1987), who noted the Seyfert- 
like characteristics of its optical spectrum, including a broad 
component to the H/3 emission line. The galaxy has a central 
velocity dispersion of 36 km s" 1 , which yields a mass esti- 
mate of Mbh = 1 -4 x 10 5 M©, again consistent with the value 
of 1.6 x 10 5 M derived from the Hf3 linewidth-luminosity- 
mass scaling relation. The [O III] A5007 linewidth, a « 37 
km s , is also fully consistent with the stellar velocity disper- 
sion (Barth et al. 2004). 

Unfortunately, it is currently technically impossible to ob- 
tain direct mass measurements for Mbh ;$ 10 6 Mq, because 
of our inability to resolve the BH sphere of influence, for 
all but the nearest galaxies in the Local Group (e.g., M33; 
Gebhardt et al. 2001). The only direct mass determina- 
tions for intermediate-mass BHs come from observations of 
massive star clusters (Gebhardt, Rich, & Ho 2002; Gerssen 
et al. 2002), which, in any case, have been controversial 
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(Baumgardt et al. 2003a, b). Barring evidence from spatially 
resolved kinematics, the best indicator for nuclear BHs in 
nearby galaxies comes from AGN activity as revealed through 
optical emission-line intensity ratios (Baldwin, Phillips, & 
Terlevich 1981; Veilleux & Osterbrock 1987; Ho, Filippenko, 
& Sargent 1993), kinematically distinct, broad emission lines 
(e.g., Ho et al. 1997c), or the presence of compact, nonstellar 
radio and X-ray cores (e.g., Ho et al. 2001; Nagar et al. 2002; 
Terashima & Wilson 2003; see Ho 2004a for a review). 

Thus, the best hope to find intermediate-mass BHs is 
through an AGN survey. However, there are at least two prac- 
tical challenges. First, it may be that objects like NGC 4395 
are quite rare. In the Palomar survey of ~ 500 nearby galax- 
ies, NGC 4395 is the only example of its kind found (Ho, 
Filippenko, & Sargent 1997a, b). POX 52 was discovered 
serendipitously in a fairly limited objective prism survey for 
emission-line objects (Kunth et al. 1987). These two are the 
only known cases of intermediate-mass BHs in galactic nu- 
clei. Thus, to accumulate a significant sample will require a 
much larger survey. Second, even radiative signatures of ac- 
cretion in these objects will be challenging to uncover. The 
host galaxies are expected to be either very late-type systems 
with low mass, and therefore intrinsically faint, or else regu- 
lar, luminous disk-dominated galaxies with a tiny or nonexis- 
tent bulge. Even at a high accretion rate a 10 5 M Q BH will 
itself have a low luminosity. As a result, detection of the faint 
AGN will be limited by sensitivity and host galaxy contami- 
nation. 

A large-area, sensitive, and uniform optical spectroscopic 
galaxy survey, such as the Sloan Digital Sky Survey (SDSS; 
York et al. 2000), offers the best opportunity for finding a sig- 
nificant number of new intermediate-mass BH candidates. We 
present the first sample of candidates culled from the SDSS. 
We discuss the data in §2, the sample definition and analysis 
in §3, and the main results in §4. In this paper we assume 
the following cosmological parameters to derive distances: 
H = 100/z = 72 km s" 1 Mpc" 1 , Q m = 0.3, and fl A = 0.7. 

2. THE SDSS DATA 

Since we have relied heavily on the SDSS software 
pipeline, we will briefly discuss relevant characteristics of the 
survey. This is by no means a complete discussion. Details 
of SDSS photometry and spectroscopy can be found in the 
Early Data Release paper (Stoughton et al. 2002; hereafter 
EDR) and the Data Release One paper (Abazajian et al. 2003; 
hereafter DR1). 

The SDSS uses a dedicated 2.5 m telescope to image one- 
quarter of the sky and obtain follow-up spectroscopy on au- 
tomatically selected star, galaxy, and quasar candidates. The 
DR1 includes a photometric sample of 53 million objects over 
2099 deg 2 and photometrically defined quasar, galaxy, and 
stellar spectra over 1360 deg 2 (Abazajian et al. 2003). 

Calibration and sample selection for both surveys is con- 
ducted by automated pipelines. The photometric pipeline, 
Photo, performs bias subtraction, flat-fielding, cosmic ray re- 
jection, and masking of bad pixels, followed by astrometry, 
point-spread function modeling, object identification, and im- 
age deblending. It measures four different (asinh) magni- 
tudes (Lupton, Gunn, & Szalay 1999). The Petrosian mag- 
nitude, measured within a radius defined by the shape of the 
azimuthally averaged light profile of the galaxy, is best for 
photometry of extended objects such as galaxies (Stoughton 
et al. 2002). Typical photometric calibration errors are 0.02 
mag in g (Smith et al. 2002), while the zeropoints between 



TABLE 1 
The SDSS Sample 
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SDSS Name 
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0.66 


0.11 


3 


*\DS^ Tfn?VlS ^04-00^408 4 

Oi^kJiJ JUJi 1 J J . J JlUO .1 
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16 
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09 
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cnss nni ins 404-007QO8 i 
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19 39 


56 


13 
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SD^S T1 01677 ^7—000714 S 

JI7JJ J 1 \J 1 UZ. / .Ji \J\J\J 1 l^t . J 


0943 


19 31 


0.66 


13 




SDSS Tl 14008 71+03071 1 4 


0.081 1 


17.32 


0.51 


0.07 


9 


SDSS J115138.24+004946.4 


0.194 


19.02 


1.20 


0.10 


10 


SDSS J124035.81-002919.4 


0.0809 


18.00 


0.47 


0.09 


11 


SDSS J125055.28-015556.6 


0.0814 


18.23 


0.37 


0.08 


12 


SDSS J135724.52+652505.8 


0.106 


18.45 


0.36 


0.07 


13 


SDSS J141234.67-003500.0 


0.126 


18.19 


0.59 


0.17 


14 


SDSS J143450.62+033842.5 


0.0281 


15.60 


0.62 


0.14 


15 


SDSS J144507.30+593649.9 


0.127 


19.98 


0.41 


0.03 


16 


SDSS J170246.09+602818.9 


0.0690 


18.22 


0.59 


0.07 


17 


SDSS J172759. 15+542147.0 


0.0994 


19.27 


0.57 


0.14 


18 


SDSS J232159.06+000738.8 


0.183 


19.39 


0.71 


0.16 


19 


SDSSJ233837.10-002810.3 


0.0355 


16.87 


0.60 


0.11 



NOTE. — Col. (1): Identification number assigned in this paper. Col. (2): Official 
SDSS name. Col. (3): Redshift measured by the SDSS pipeline. Col. (4): Petrosian g 
magnitude. Col. (5): Petrosian g - r color. Col. (6): Galactic extinction in the g band. 



the SDSS and AB magnitude system could differ by as much 
as 10% (Stoughton et al. 2002). 

Galaxy targets are automatically selected based on mor- 
phology. Specifically, resolved sources with r < 17 '.77 mag 
are targeted as galaxies (note that this cut-off varied between 
the EDR and the DR1). Candidates with magnitudes brighter 
than 15 in g or r or 14.5 in i are rejected in case of contam- 
ination of other fibers. A surface brightness cut is also ap- 
plied, and targets may be excluded due to proximity because 
the fiber holders have a finite size and are kept 55" apart from 
each other to avoid fiber collisions. Quasar candidates are se- 
lected on the basis of a color selection technique detailed in 
Richards et al. (2002). 

SDSS spectra are acquired with a pair of double, fiber-fed 
spectrographs. A plate of 3 deg 2 with 640 optical fibers is 
drilled for each field. Each fiber subtends a diameter of 3", 
corresponding to ~ 6.5 kpc at z = 0.1. The spectra have an 
instrumental resolution of A/AA s=a 1800 (cr ms « 71 km s" 1 ). 
Integration times are determined for a minimum signal-to- 
noise ratio of 4 at g = 20.2 mag. The spectroscopic pipeline 
performs basic calibration as above, as well as spectral ex- 
traction, sky subtraction, removal of atmospheric absorption 
bands, and wavelength and spectrophotometric calibration. 
Redshifts are determined using x 2 fits to stellar, emission-line 
galaxy, and quasar templates, as well as to the emission lines 
themselves; the typical accuracy is ~ 30 km s" 1 . 
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FIG. 1. — Left: Spectrum with (bottom) and without (top) starlight and continuum subtraction. For clarity, the top spectrum in some cases is offset by 
an additive constant, which is indicated. The top left corner of each panel gives the object identification number, which can be referenced in Table 1. Right: 
Enlargement of the region around Ha+[N n] AA6548, 6583. In the upper plot, the data are shown in histogram, the individual models for the narrow and broad 
lines are shown in dotted lines, and the sum of all the fitted components is shown as a solid line. The lower plot gives the residuals of the fit. The ordinate of the 
plots are in units of 10~ 17 erg s -1 cm -2 A -1 . 
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Fig. 2.— 



Same as Figure 1. 
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Same as Figure 1 
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3. THE SAMPLE AND ANALYSIS 

Our primary sample was culled from a parent sample of 
broad-line AGNs at z < 0.35 in the DR1. These are un- 
ambiguous AGN candidates with virial BH mass estimates 
< 10 6 Mq, which hereafter we consider to be the upper bound 
for intermediate-mass BHs. We determine "virial" BH masses 
using a variant of the linewidth-luminosity-mass scaling rela- 
tion of Kaspi et al. (2000), whose robustness has been tested 
with respect to theM B H-cr* relation of inactive galaxies (Geb- 
hardt et al. 2000b; Ferrarese et al. 2001; see Barth 2004 for a 
review). In order to apply this method, we first identify all the 
AGNs from DR1 that have a broad component to at least the 
Ha emission line. In order to detect objects with weak lines 
and to reliably measure emission-line and continuum param- 
eters, we must properly remove the contaminating starlight 
from the host galaxy, model the profiles of the emission lines, 
and quantify the strength of the AGN continuum luminosity. 
We discuss these various steps below. 

3.1. Starlight and Continuum Subtraction 

Low-luminosity AGNs are challenging to detect because 
the signature emission lines are often contaminated by 
starlight. The SDSS fibers always include significant host 
galaxy emission, making careful removal of the stellar contin- 
uum essential to reliably measure the emission-line spectrum, 
and especially to detect any weak, broad features. Because 
we expect Ha to be the strongest permitted line in the optical, 
and hence best suited for detecting a broad-line component if 
present, we start with the <~ 153,000 galaxies and quasars in 
the DR1 with z < 0.35. 

The next step is to remove the starlight, along with any ad- 
ditional featureless nonstellar continuum emission, in order to 
produce a pure emission-line spectrum. A variety of template- 
fitting techniques exist to accomplish this (see, e.g., Ho 2004a 
for a review). For this work we use a principal component 
analysis (PCA) algorithm adopted from Hao & Strauss (2004; 
see also Connolly & Szalay 1999). Assuming a linear rela- 
tion between all spectral variations, PCA (or the Karhunen- 
Loeve transform) transforms a large library of input spectra 
to an orthogonal basis in which each eigenspectrum is along 
a direction of maximal variance. The input spectra are de- 
rived from absorption-line spectra from the SDSS, and an A- 
star spectrum may be included to account for a post-starburst 
stellar population. Each program spectrum is then modeled 
as a linear combination of the eigenspectra. The PCA tech- 
nique has a number of benefits when a large, homogeneous 
database like SDSS is available. The method typically re- 
quires no more than eight eigenspectra to fully reconstruct a 
given galaxy spectrum, and for absorption-line galaxies the 
resulting solution is unique. Furthermore, the eigenvectors 
naturally include variance in velocity dispersion as long as 
the template galaxies span the same range as the target sam- 
ple. In practice, some galaxies contain substantial nonstellar 
continuum emission, which we model as a power law fixed 
to a slope of (i = 1.0, where f\ oc \~P . As discussed in J. 
E. Greene & L. C. Ho (in preparation), the AGN-dominated 
objects in our sample typically have a power-law continuum 
spectrum with this slope. 

Examples of spectra before and after starlight and contin- 
uum subtraction are given in the left panels of Figures 1-4. 

3.2. Initial Selection of Broad Ha Candidates 

Once we obtained continuum-subtracted spectra, we used, 
as a first pass, a simple algorithm to cull objects with ex- 



cess emission in the region 6300-6700 A that is potentially 
broad Ha emission. We began by removing any narrow emis- 
sion lines from this spectral region by adaptively smoothing 
it, iteratively removing outlier points. We then calculated the 
mean, root mean square (rms) normalized to the dispersion in 
a line-free spectral region (5600-6000 A), and integrated flux 
of the region. Since the spectra are continuum-subtracted, any 
object with a nonzero mean, rms, or integrated flux is a can- 
didate for having broad Ha emission. By trial and error we 
chose thresholds in the mean, integrated flux, and normalized 
rms that identified all broad-line objects down to a very low 
contrast. We conservatively chose a low threshold, which in- 
creased our contamination rate but enhanced our complete- 
ness. The resulting sample from this procedure comprised 
4163 candidates. As an initial cut, this technique has the ben- 
efit of being computationally inexpensive (important for an 
initial sample as large as the DR1). It has two major disad- 
vantages: potential incompleteness and high contamination. 
Contaminants included misclassified late-type stars, which 
tend to have a broad maximum near 6500 A that can be mis- 
taken for broad Ha emission (Ho et al. 1997c), and objects of 
unknown origin with large calibration errors or whose spec- 
tra have been corrupted in the Ha region. We experimented 
with using pipeline flags to automatically remove these, but 
in the end flagged them by hand. A more significant con- 
taminant comes from incomplete removal of narrow emission 
lines from starburst galaxies with extremely strong line emis- 
sion, or emission-line galaxies with multi-component, rela- 
tively broad narrow lines. Detailed profile fitting, as described 
below, is required to isolate galaxies with a truly kinemati- 
cally distinct broad-line component. 

3.3. Emission-line Profile Fitting 

The broad Ha feature, a reliable signature of AGN activ- 
ity (see Ho et al. 1997c), is coincident with a combination 
of narrow Ha emission and the flanking [N II] AA6548, 6583 
doublet. While in many cases the broad-line emission in our 
candidates is evident by eye, the Ha+[N II] region is suffi- 
ciently complicated that some care is needed to determine the 
presence or absence of a broad Ha component. In addition to 
imperfect removal of starlight, the shape of the narrow lines 
themselves may confuse the detection of broad Ha. There are 
often broad, non-Gaussian wings, multiple peaks, and asym- 
metries in the narrow lines that can mimic or hide a true broad 
component to Ha. We therefore use a profile-fitting proce- 
dure, similar to that of Ho et al. (1997c), to model the narrow 
components based on regions of the spectrum less prone to 
blending and confusion. We can then determine objectively 
whether a broad component is present. 

While the profiles of the narrow emission lines in some 
AGNs are known to vary as a function of critical density and 
ionization parameter (e.g., Pelat, Alloin, & Fosbury 1981; 
Filippenko & Halpern 1984; Filippenko 1985; De Rober- 
tis & Osterbrock 1986), empirically the line profile of [S II] 
AA6716, 6731 is generally well matched to those of [N II] and 
narrow Ha (Filippenko & Sargent 1988; Ho et al. 1997c). 
Therefore, we model the [S II] doublet with a superposition 
of Gaussian components and then scale this model to fit the 
Ha+[N II] narrow lines. We use up to four Gaussian com- 
ponents when it is statistically justified. We do not ascribe 
any particular physical significance to the individual compo- 
nents, but they allow us to describe the line shape as fully 
as the data quality allows. The two [S II] lines are assumed 
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FIG. 4. — Same as Figure 1. 

to have the same profile and are fixed in separation by their 
laboratory wavelengths; the relative strengths of the two lines 
change with density and so are allowed to vary in the fit. In a 
few cases where [S II] is undetectable or too weak to yield a 
reliable fit, we construct a model using [O III] A5007. [O III] 
is both strong and nearly ubiquitous, and may seem like a bet- 
ter choice for template fitting. Unlike [S II], the [O III] profile 
does not typically match that of [N II] in detail. [O III] of- 
ten has a broad, blue shoulder to its profile, suggestive of an 
outflow origin (De Robertis & Osterbrock 1984, 1986; Whit- 
tle 1985). Nonetheless, [O III] still seems to be an acceptable 
substitute for [S II] when the latter is unavailable. 

Once a satisfactory narrow-line model is achieved it serves 
as a template for fitting the Ha+[N II] complex. The centroids 
of each component are fixed relative to each other at their lab- 
oratory separations, and the relative strength of [N II] A6583 
and [N II] A6548 is fixed at the theoretical value of 2.96. An 
additional broad component for Ha is added if statistically 
warranted; we use as many Gaussians as needed, but gener- 
ally no more than 2-3, to fully describe the complex shapes of 
some of the broad components. Our fits are shown as dotted 
lines in the right panels of Figures 1-4. From the final model 
for the broad Ha line, we measure the full width at half maxi- 
mum (FWHM) of the profile and its integrated flux. This step 
eliminated a large number of broad Ha candidates, leaving a 
sample of ^3200 genuinely broad-line objects. 

3.4. Nonstellar Continuum Luminosity 

As discussed below, in order to estimate BH masses we 
need to know the luminosity of the featureless nonstellar 
(AGN) continuum, L5100 = ALa at A = 5100 A. We obtain 
L5100 in two ways. For objects whose AGN component clearly 
dominates (e.g., object #1, 5, and 6 in Fig. 1), we directly 
fit a single power law to the spectrum using spectral regions 
known to be relatively unaffected by strong emission lines, 
especially from Fe II, which is pervasive in AGN spectra and 
can often be mistaken for continuum emission. The "contin- 
uum" windows we chose, guided by Francis et al. (1991), 
are: 3700-3715, 3740-3800, 4041^tt)43, 4150-4250, 5550- 
5850, 6000-6290, and 6400-6450 A. From the power-law fit 
we derive a slope and the continuum flux at 5100 A. 

For spectra with noticeable galaxy contribution, we must 
rely on the PCA decomposition to model the power-law con- 
tribution. From extensive experimentation with the PCA fits, 
we find that we can reliably recover the power-law component 
in objects where the AGN contribution at 5100 A is > 20%. 
In this study we only consider objects whose L5100 exceeds 
this threshold. 

Once we have a measurement of the AGN power-law con- 
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tinuum, we can translate it into an equivalent g-band abso- 
lute magnitude, M^(AGN), by convolving the power law with 
the response function of the SDSS g filter. In turn, we can 
estimate the g-band absolute magnitude of the host galaxy, 
M^(host), by subtracting the AGN contribution to the total ob- 
served g-band luminosity. 

3.5. Black Hole Masses 

We estimate BH "virial" masses using a variant of 
the linewidth-luminosity-mass scaling relation discussed by 
Kaspi et al. (2000). If the broad-line region (BLR) is in 
virial equilibrium, then the mass of the central BH is M B h~ 
w 2 /?blr/G, where v = (V3/2) x wfwhm and wfwhm is the 
FWHM of the broad component of H/3, Rblr is the radius of 
the BLR, and G is the gravitational constant. For a sample of 
34 low-redshift type 1 AGNs studied with reverberation map- 
ping Kaspi et al. (2000) found an empirical relation between 
the continuum luminosity of the AGN at 5100 A and the size 
of the BLR: 



Ablr = (32.9!^) 



^5100 



lO^ergs" 1 

With this scaling relation, the virial mass is simply 



0.700±0.033 



lt-days. (1) 



M BH = 4.82x 10 6 



Lsi 



10 44 erg s" 



/ WFWHM 

V 1000km s 



Mr,. 



(2) 

While it is conventional to use H/3 to determine vfwhm, we 
have decided to use Ha as a surrogate instead. Ha is at least 
a factor of 3 stronger than H/3, and therefore much easier to 
measure for low-luminosity sources or in objects where host 
galaxy contamination is severe. Another advantage of Ha is 
that, despite blending with [N II] emission, it is still much 
easier to measure because it is not strongly affected by nearby 
Fe II emission, as is the case for H/3 (see, e.g., Francis et al. 
1991). 

The currently popular formalism to derive virial BH masses 
for broad-line AGNs is calibrated against the sample of AGNs 
with reverberation mapping data (Kaspi et al. 2000). As 
discussed by Krolik (2001), this simplistic approach may be 
subject to a variety of unquantified systematic uncertainties. 
Nevertheless, Gebhardt et al. (2000b) and Ferrarese et al. 
(2001) have demonstrated, for a subset of the reverberation- 
mapped AGNs with available stellar velocity dispersions, that 
the virial BH masses, remarkably if unexpectedly, do roughly 
conform to the M Bli ~ 0* relation of inactive galaxies. Interest- 
ingly, both NGC 4395 and POX 52 follow the low-mass ex- 
trapolation of the Tremaine et al. (2002) fit to the M B h-o* re- 
lation (Filippenko & Ho 2003; Barth et al. 2004; Ho 2004b). 
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TABLE 2 
Emission-line Measurements 



ID 


[On] 


Fe ii 


(H/3)„ 


(H/3) b 


[O in] 


[Oi] 


(Ha)„ 


(H«) b 


[N n] 


[Sn] 


[S n] 


rWHM Ha 


FWHM [0 m] 




A3727 


A4570 






A5007 


A6300 






A6583 


A6716 


A6731 


(km s ') 


SI I N 

(kms ') 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


1 


154 


1 SQ 


046 


1 48 


362 


021 


143 


A 73 


n HQS 






R30 




2 


273 


<1.04 


079 


88 


519 


053 


245 


1 95 


0.141 


077 


074 


732 


262 


3 


0.225 


1 07 


0.070 


0.73 


192 


<0 002 


217 


3 25 


051 


064 


038 


916 


199 


4 


0.307 


0.26 


139 


0.42 


1880 


030 


431 


1 48 


0.165 


0.073 


063 


870 


238 


5 


081 


4 29 


167 


1.79 


210 


<0 005 


517 


7.98 


0.400 


0.097 


0.055 


698 


291 




171 


2 83 


074 


1 Q3 


104 


<0 008 


461 


8 31 


206 


099 


u.uo / 


Q1 Q 

717 






309 


55 


093 


25 


453 


08 1 


288 


1 16 


170 


090 


095 


940 


449 




0.096 


6.98 


0.311 


1.89 


241. 


0.182 


0.966 


11.0 


1.00 


0.242 


0.241 


591 


297 


9 


0.085 


<0.55 


0.019 


0.07 


1050. 


<0.001 


0.060 


0.46 


0.062 


0.014 


0.017 


703 


504 


10 


0.240 


0.41 


0.127 


0.25 


2150. 


0.048 


0.396 


1.35 


0.138 


0.071 


0.069 


722 


255 


11 


0.578 


0.44 


0.219 


0.20 


1140. 


0.043 


0.680 


1.01 


0.203 


0.155 


0.118 


642 


212 


12 


0.133 


0.92 


0.154 


0.75 


463. 


0.020 


0.478 


2.87 


0.263 


0.040 


0.054 


872 


231 


13 


0.191 


3.55 


0.157 


1.58 


126. 


<0.007 


0.487 


6.65 


0.335 


0.114 


0.074 


785 


275 


14 


0.219 


<1.07 


0.273 


0.52 


541. 


0.053 


0.848 


1.96 


0.523 


0.159 


0.140 


1089 


276 


15 


0.079 


2.12 


0.012 


2.93 


83.4 


<0.()04 


<0.()38 


11.8 


<0.038 


<0.061 


<0.061 


750 


219 


16 


0.565 


<1.53 


0.211 


1.18 


257. 


0.028 


0.654 


6.17 


0.124 


0.148 


0.115 


919 


234 


17 


0.244 


0.97 


0.028 


0.65 


386. 


0.033 


0.086 


2.08 


0.146 


0.059 


0.051 


464 


233 


18 


0.707 


<7.30 


0.325 


0.75 


49.5 


<0.038 


1.01 


7.24 


0.899 


<0.434 


<0.434 


707 


342 


19 


0.234 


<0.13 


0.075 




2550. 


0.045 


0.232 


0.16 


0.175 


0.072 


0.069 


1730 


216 



NOTE. — Col. (1): Identification number assigned in this paper. Col. (2)-(12): All fluxes are relative to that of [O ill] A5007, which is listed in units of l(T 17 erg s~' cm" 2 . Note 
that these are observed values; no extinction correction has been applied. The subscripts "n" and "b" in Col. (4)-(5) and (8)-(9) refer to the narrow and broad components of the line, 
respectively. The narrow component of H/3 is sufficiently weak that it was impossible to obtain a robust measurement of its strength, so we fixed it to that of the narrow Ha component, 
assuming an intrinsic ratio of Ha/H/3 = 3.1 for AGNs (Halpern & Steiner 1983; Gaskell & Ferland 1984). Col. (13)— (14): These are the observed linewidths; they have not been 
corrected for instrumental resolution. The Ha FWHM was derived from the model fit to the broad Ha line, as described in §3.3. 



The virial masses scatter around the M B h - o* relation with a 
spread of ~ 0.5 dex, which for the moment can be taken as an 
estimate of the systematic uncertainty of this method. 

Certain aspects of our analysis may introduce additional er- 
rors into the mass measurements, although we believe that 
they are small relative, or at most comparable, to the intrin- 
sic systematic uncertainty discussed above. There are two ba- 
sic measurements that affect our mass determinations, namely 
L5100 and ufwhm- Both of these quantities become increas- 
ingly challenging to measure when galaxy contamination is 
large. The complete profile of weak emission lines with very 
broad wings becomes very difficult to measure when the con- 
trast between the AGN and the galaxy is low. Dilution from 
host galaxy starlight, which can never be removed perfectly, 
always leads one to underestimate the true width of a weak, 
low-contrast broad emission line. This has been demonstrated 
most conclusively for some objects by comparing starlight- 
subtracted ground-based spectra with small-aperture spectra 
taken with the Hubble Space Telescope (e.g., Ho et al. 2000; 
Shields et al. 2000): the small-aperture spectra, largely free 
from starlight contamination, reveal much broader Ha wings 
than the ground-based spectra. To mitigate this problem, and 
that associated with measuring L5100 (§3.4), in this study we 
have chosen to concentrate on objects where the AGN com- 
ponent is largely dominant. This is evident from Figures 1^4- 
by the overall weakness of stellar absorption features in the 
spectra. 

We note that substituting Ha for H/3 to measure ufwhm 
should introduce negligible additional uncertainty into the 
mass estimates. Although the physical conditions in the BLR 
that give rise to Ha and H/3 emission in principle may be 
slightly different (e.g., Kwan & Krolik 1981), in practice the 
two lines are observed to have rather similar velocity pro- 
files, with H/3 on average being only slightly broader than Ha. 
From a study of 19 Seyfert 1 galaxies, Osterbrock & Shuder 



(1982) find that the broad H/3 and Ha lines have nearly iden- 
tical full widths at zero intensity, and on average H/3 is 16% 
broader than Ha at FWHM. If this trend holds, then we may 
be underestimating the mass in Equation 2 by ~30%, which is 
insignificant compared to the much larger 0.5 dex uncertainty 
in the virial mass zeropoint. Moreover, in the two AGNs with 
intermediate-mass BHs known to date — NGC 4395 and POX 
52 — there is no indication that the line profile of broad Ha 
differs systematically from that of H/3 (Filippenko & Sargent 
1989; Kraemer et al. 1999; Barth et al. 2004). We await future 
measurements of stellar velocity dispersions to investigate po- 
tential systematic differences in our virial mass estimates, but 
for now we assume that the Ha linewidth can be substituted 
for that of H/3. 

With the BH virial masses at hand, we are in a position to 
select the intermediate-mass BH candidates. With a cut-off of 
10 6 Mq, the final sample contains 19 broad-line objects (Fig. 
5). Their spectra are in Figures 1-4. Table 1 shows the basic 
properties of the sample, Table 2 lists measurements of the 
emission-line parameters, and Table 3 gives the luminosity 
and mass estimates. 

4. DISCUSSION 

We employ a simple technique to exploit the breadth of the 
SDSS galaxy spectroscopy to extend the demography of cen- 
tral BHs at least 1 order of magnitude below the 10 6 M Q 
threshold, a regime hardly explored previously. We present 
the first sizable sample of AGNs with candidate intermediate- 
mass BHs. The sample contains 19 objects with virial BH 
masses ranging fromM B H ~ 8 x 10 4 to 10 6 M Q (Fig. 5). We 
now briefly examine some of the properties of these objects. 

4. 1 . Spectral Properties 

In the literature to date, the closest analog to our sample 
(aside from NGC 4395 and POX 52) is the class of AGNs 
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TABLE 3 

Luminosity and Mass Measurements 



ID 


M g (total) 


M g (AGN) 


M g (host) 




logL 510 o 


log M BH 


logL bol /L Edd 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


1 


—19.1 


—17 8 


—18 7 


60 


43.06 


5 86 


092 


2 


_i o q 
lo.y 


— 1 S 4 


_1 Q Q 

lo.y 


1 on 


HI .o / 


4 92 


—A 1 S 

KJ. 13 


3 


—18.7 


—17.4 


—18.4 


0.11 


42.91 


5.84 


—0.03 


4 


—18.6 


—17 1 


—18.3 


0.88 


42.80 


5.72 


—0.02 


5 


—19.2 


—17 7 


—18.9 


0.27 


43.10 


5.74 


0.254 




—18.8 


—17.7 


—18 3 


0.076 


42.93 


5 86 


—0.03 


7 


—18.8 


—17 1 


—18.5 


0.44 


42.97 


5.91 


—0.04 


8 


-20.4 


-18.3 


-20.2 


0.14 


43.35 


5.77 


0.475 


9 


-20.8 


-18.4 


-20.6 


1.00 


43.02 


5.69 


0.223 


10 


-19.7 


-18.4 


-19.4 


0.46 


43.33 


5.93 


0.296 


11 


-19.5 


-17.5 


-19.3 


0.76 


43.22 


5.75 


0.362 


12 


-19.9 


-18.1 


-19.7 


0.79 


43.17 


5.98 


0.082 


13 


-20.6 


-18.1 


-20.5 


0.39 


43.32 


5.99 


0.218 


14 


-19.7 


-15.4 


-19.7 


1.00 


41.92 


5.30 


-0.48 


15 


-18.8 


-18.3 


-17.8 


0.24 


42.97 


5.71 


0.153 


16 


-19.1 


-17.4 


-18.9 


1.00 


42.53 


5.58 


-0.15 


17 


-18.9 


-17.6 


-18.5 


0.31 


43.05 


5.35 


0.592 


18 


-20.3 


-18.0 


-20.1 


1.00 


43.09 


5.74 


0.240 


19 


-19.0 


-15.0 


-19.0 


1.00 


41.69 


5.54 


-0.95 



NOTE. — Col. (1): Identification number assigned in this paper. Col. (2): Total g-band absolute magnitude. Col. (3): AGN g-band absolute magnitude, estimated from L5100 given 
in Col. (6) and a conversion from L5100 to M s assuming f\ oc A - * 9 , where /3 is given in Col. (5); typical uncertainty is ~ 0.4 mag. Col. (4): Host galaxy g-band absolute magnitude, 
obtained by subtracting the AGN luminosity from the total luminosity. Col. (5): Power-law slope /3, where f\ oc A - * 9 ; typical uncertainty is ~ 10%. A value of 1.00 indicates that it 
was fixed in the PCA decomposition. Col. (6): AGN continuum luminosity at 5100 A (erg s~'), estimated from a power-law fit to the continuum; typical uncertainty is ~ 10%. Col. 
(7): Virial mass estimate of the BH (Mq). Col. (8): Ratio of the bolometric luminosity (see text) to the Eddington luminosity. 



ranging from 464 to 1730 km s , with an average value of 
836 km s -1 . However, we note that our sample shows a greater 
diversity of Fe II and [O III] strengths than in previous sam- 
ples of NLS Is. Fe II emission is common, but not ubiquitous 1 . 
For example, the intensity ratio of Fe II A4570 relative to total 
(broad + narrow) H/3, where Fe II A4570 is the complex be- 
tween 4434 and 4680 A, varies from ~0.465 to 3.17, with 
a mean of 1.32 ±0.16, calculated using the Kaplan-Meier 
product-limit estimator to account for upper limits (Feigel- 
son & Nelson 1985). This is to be compared with (Fe II 
A4570/H/3(total) ) = 0.67 ± 0.04 for the 56 NLS Is presented 
by Veron-Cetty, Veron, & Goncalves (2001). Similarly, in 
our sample [O III]/H/3(total) = 0.335 to 12.0, with a mean 
of 2.24 ±0.72. In contrast Veron-Cetty et al. (2001) find 
([O III]/H/3(total)) = 0.65 ±0.09, and Constantin & Shields 
2003 find ([O III]/H/3(total)) = 0.39 for their composite spec- 
trum of 22 NLS Is. Thus, compared to classical NLS Is, we 
find stronger [O III] and Fe II emission. The distribution of 
power-law slopes for the AGN-dominated spectra are also 
quite interesting (Table 3). While 1.5 < (i < 1 is a standard 
range for luminous AGNs (e.g., Malkan 1988), we find an av- 
erage value of (3 = 0.42. 

Only six of the 19 sources (30%) were detected in the 
ROSAT All-Sky Survey (Table 4), confirming the finding of 
Williams et al. (2002) that not all NLS Is have detectable 
soft X-ray (0.5-2 keV) emission. The data are not of suffi- 
ciently high quality to make spectral analysis possible, so to 
convert count rate to flux we assume a power-law spectrum 
with a fixed photon index of T = 3, which is roughly the av- 
erage value observed in previous ROSAT studies of NLS Is 

1 The same comment applies to NGC 4395 and POX 52, both of which 
may be formally regarded as NLSls. Fe ii multiplets can be seen in the UV 
spectrum of NGC 4395 (Filippenko et al. 1993), but they are quite weak in 
its optical spectrum (Filippenko & Sargent 1989). Fe n emission is also not 
prominent in the optical spectrum of POX 52 (Barth et al. 2004). 
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FIG. 5. — The distribution of virial BH masses for our sample. 



known as narrow-line Seyfert Is (NLS1), which are thought 
to possess relatively low-mass BHs emitting at a high fraction 
of their Eddington rates (e.g., Boroson 2002). Since our sam- 
ple selection required a relatively prominent AGN component 
relative to the host galaxy, the sample should preferentially 
single out AGNs in a relatively "high" state. 

First classified as a group by Osterbrock & Pogge (1985), 
NLSls are typically characterized by broad permitted lines 
with FWHM < 2000 km s" 1 , prominent Fe II emission, and 
relatively weak [O III] AA4959, 5007 lines. NLSls also tend 
to emit copiously in the soft X-ray band (e.g., Boiler, Brandt, 
& Fink 1996; but see Williams, Pogge, & Mathur 2002), 
presumably as a consequence of their high accretion rates 
(Pounds, Done, & Osborne 1995). Our sample conforms to 
some, but not all, of these expectations. By the linewidth cri- 
terion, certainly all of the objects would technically qualify as 
"narrow-line" sources: the broad Ha component has FWHM 
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FIG. 6. — Diagnostic diagram plotting log [O in] A5007/H/3 versus log [N ii] A6583/Ha (a) and versus log [S n] AA6716, 6731/Ha (b). The line ratios have 
not been corrected for reddening, but this should not matter because of the close wavelength separation of the lines. Object #15 was omitted because [N II], [S 
II], and narrow Ho were not detected. The dotted lines mark the boundaries of the three main classes of emission-line nuclei, according to the convention of Ho 
et al. (1997a): "H" = H ii nuclei, "S" = Seyferts, and "L" = LINERs. 



TABLE 4 
X-ray Detections 



ID 


C 


log Ah 


log/x 


logZ-x 


(1) 


(2) 


(3) 


(4) 


(5) 


1 


0.0483 


20.59 


-12.21 


42.92 


7 


0.0276 


20.59 


-12.46 


42.86 


8 


0.113 


20.37 


-12.86 


42.32 


12 


0.100 


20.30 


-12.92 


42.51 


14 


0.0564 


20.40 


-12.16 


42.07 


17 


0.0174 


20.50 


-12.66 


42.70 



NOTE. — Col. (1): Identification number assigned in this paper. Col. (2): ROSAT 
count rate (counts s _1 ). Col. (3): Galactic column Nh (cm""), calculated following 
Dickey & Lockman (1990). Col. (4): X-ray flux in the 0.5-2 keV band (erg s~ ] m 2 ) 
assuming a power-law spectrum with T = 3 and Nh from Col. (3). Col. (5): X-ray 
luminosity in the 0.5-2 keV band (erg s~ ). 



Figure 6 places the sources on two commonly used line in- 
tensity ratio diagnostic diagrams. The dotted lines schemat- 
ically denote the boundaries used by Ho et al. (1997a) to 
demarcate the main classes of emission-line objects. Al- 
though all the objects are clearly AGNs by the presence of 
their prominent broad lines, some of them in the [O III]/H/3 
vs. [N H]/Ha and [O HI]/H/3 vs. [S Il]fHa diagrams formally 
lie outside of the AGN zone ([N II]/Ha > 0.6 and [S U]fHa 
> 0.4). This is rarely seen in local type 1 AGNs (Ho et al. 
1997a). Since the SDSS fibers subtend a significant fraction 
of each galaxy, the narrow emission-line ratios of the nuclei 
are, at least in part, contaminated by host galaxy emission, 
which would have the effect of lowering the [S II]/Ho! and 
[N II]/Ha ratios. Thus, one expects, a priori, that AGNs not 
selected by their narrow lines will, in some cases, exhibit 
narrow-line ratios that fall below the AGN threshold. 



(e.g., Boiler et al. 1996; Veron-Cetty et al. 2001; Williams 
et al. 2002). We use the online program webPIMMS 2 to 
do the conversion, assuming that the only source of absorp- 
tion is due to the Galaxy along the line of sight. The de- 
rived soft X-ray luminosities range from ~ 10 42 1 to 10 42 9 
erg s" 1 (Table 4), which probably can be attributed mostly to 
the AGN. The host galaxies of our sample are relatively faint 
(§4.3), with absolute g-band magnitudes that roughly trans- 
late to L B ps 10 9 - 4 - 10 10 X L Q (Fukugita et al. 1996). If these 
galaxies follow the L B -L X relation of normal galaxies (e.g., 
Fabbiano 1989), they are expected to have L x ~ 10 38 5 erg s" 1 . 

Only one of the objects was detected in the VLA FIRST 
survey (Becker, White, & Helfand 1995), although the entire 
sample overlapped with the FIRST survey area. Object #10 
has a 20 cm flux density of 1.3 mJy, which corresponds to 
a spectral power of 2.0 x 10 22 W Hz" 1 . This level of radio 
emission is not atypical of normal galaxies (Condon 1992), 
so we cannot conclude, without additional evidence, that it is 
associated with AGN activity. 

2 |http : / /heasarc ■ gsf c ■ nasa ■ gov/ Tool s/w3pimms . html| 



4.2. Nuclear Luminosities and Eddington Ratios 

Since our sample was selected to have low BH mass, we 
expect that the AGNs will generally be of relatively low lu- 
minosity. Figure 7 shows the distribution of nuclear g-band 
absolute magnitudes for the sample. To estimate this quantity, 
we first compute L5100 (as described above), which is then 
transformed to M g . For those spectra with direct power-law 
fits, we use the best-fit spectral index, while for the remain- 
der we assume an f\ oc A -1 spectrum. The AGN component 
has absolute magnitudes ranging fromMg pa -15 to -18 mag. 
Note that these nuclear luminosities, though generally not as 
low as many local Seyferts (Ho & Peng 2001), are still sig- 
nificantly lower than in many previously studied low-redshift 
type 1 AGNs (e.g., Kohleret al. 1997). 

It is of interest to compare the nuclear luminosities with re- 
spect to the Eddington luminosities. The selection criteria — 
low BH mass and a prominent AGN spectrum — are expected 
to favor objects with high Eddington ratios. A major uncer- 
tainty in this exercise is that we do not have a robust estimate 
of the bolometric luminosity, but only a measurement of the 
optical luminosity (Lgioo)- Since the spectral energy distribu- 
tions of AGNs are known to vary significantly with accretion 
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FIG. 7. — Distributions of g-band absolute magnitudes. The magni- 
tudes pertain to the entire system (top), the host galaxy alone (middle), and 
the AGN alone (bottom). The AGN contribution was estimated from L 5m , 
which was then converted to a g magnitude assuming an f\ oc \~P spectrum, 
where is shown in Table 3. 



rate (Pounds et al. 1995; Ho 1999), it is doubtful that a single 
bolometric correction is always applicable. Nonetheless, as- 
suming Lboi = 9.8L5100 (McLure & Dunlop 2004), the derived 
Eddington ratios mostly cluster around Lboi/^Edd ~ 1 (Fig. 8; 
Table 3), as expected. 

4.3. Properties of the Host Galaxies 

The morphology of the host galaxies of intermediate-mass 
BHs is of fundamental importance for understanding the ori- 
gin of this class of objects and their relationship to the overall 
demography of central BHs in galaxies. The two currently 
known nuclear intermediate-mass BHs are found in anoma- 
lous hosts compared to the local AGN population (Ho et al. 
1997b, 2003; Kauffmann et al. 2003), which is invariably 
affiliated with massive, bulge-containing galaxies. Similarly, 
essentially every galaxy known to contain a supermassive BH 
has a bulge. Yet, NGC 4395, which without a doubt lacks a 
bulge, evidently hosts a nuclear BH (Filippenko & Ho 2003). 
The same could be said of POX 52 (Ho 2004b). Although 
POX 52 appears to be a dwarf elliptical galaxy (Barth et al. 
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FIG. 8. — Distribution of Eddington ratios. The bolometric luminosity was 
estimated from L5100 (see text). 



2004), this class of spheroids should be regarded as physi- 
cally distinct from classical bulges (Bender, Burstein, & Faber 
1992), and they may have originated as tidally stripped late- 
type disk galaxies (e.g., Moore et al. 1996). Thus, both NGC 
4395 and POX 52 provide compelling evidence that a classi- 
cal bulge is not necessary for the formation of a central BH. 

We have attempted to estimate the luminosities of the host 
galaxies by subtracting the contribution of the AGN from the 
total (Petrosian) g-band magnitude of the AGN plus host. As 
illustrated in Figure 7 (see also Table 3), in most cases the lu- 
minosities of the AGNs are sufficiently low that they hardly 
perturb the distribution of total magnitudes; thus, the host 
galaxy magnitudes should be quite robust. Interestingly, the 
hosts are relatively low-luminosity and therefore more likely 
to be late-type galaxies: Mg(host) ranges from -17.8 to -20.6 
mag, with an average value of -19.2. According to Blanton 
et al. (2003), the luminosity function of galaxies at z = 0.1 
has M* = -19.39 + 51og/z. For h = 0.72, as assumed here, 
M* = -20.1 mag. Thus, the hosts of our sample are ~ 1 mag 
fainter than M*. 

Our selection technique is unbiased with respect to galaxy 
morphology. Unfortunately, we have very limited informa- 
tion on the morphological types of the host galaxies. The low 
resolution of the SDSS imaging makes morphological clas- 
sification difficult for these objects. By selection, the galaxy 
centers are dominated by an AGN core, which makes it im- 
possible to say anything meaningful about the central struc- 
ture of the galaxies. If the low BH masses are accompanied 
by a low-mass (low-luminosity) galaxy or bulge, this makes 
it doubly hard to detect. We show a selection of the more in- 
teresting images in Figure 9. Based on the SDSS imaging, 
one-third of the galaxies clearly have a disklike component. 
Of the six that show disks, five have significant galaxy light in 
their spectra, which required PCA decomposition. Given the 
low luminosities of these galaxies, perhaps they are analogs 
of NGC 4395. The remainder are barely resolved enough to 
tell, but could be consistent with compact spheroids. We cau- 
tion that the images are of low resolution and sensitivity, so in 
many cases we may be insensitive to an extended stellar com- 
ponent. Higher resolution, deeper imaging, for example with 
the Hubble Space Telescope, is essential to determine the true 
morphology and detailed structural parameters of these hosts. 

In lieu of resolved morphological information, the SDSS 
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FIG. 9. — Sample SDSS g-band images of six objects. The vast majority 
of the sample are barely resolved and look similar to object #1. In all images 
North is up and East is to the left. Each panel is 30" X 30". 



collaboration has developed statistical proxies that describe 
morphological trends (Shimasaku et al. 2001; Strateva et 
al. 2001). We consider the "eclass" (a PCA statistic), the 
u — r color, and the "inverse concentration index." The eclass 
ranges from -0.35 for early-type galaxies to +0.5 for late-type 
galaxies. The mean eclass of our sample is 0.24. This statis- 
tic is difficult to interpret because of the strong contamination 
by the AGN. For ~ 150,000 SDSS galaxies, Strateva et al. 
(2001) find a bimodal distribution mu-r color, with the di- 
vision at u — r = 2.22 mag, which roughly divides early and 
late types. Our galaxies have a mean color of u — rK> 1.51, 
consistent with early-type galaxies. An inverse concentration 
index of ~ 0.3 corresponds to a de Vaucouleurs (1948) r 1 / 4 
profile, while an exponential profile has an index of ~ 0.43. 
The mean inverse concentration index of our objects is 0.41. 

4.4. Implications and Future Work 

While the discovery of the M BH - cr* relation established a 
compelling link between the evolution of galaxies and cen- 
tral BHs, models of this process are faced with significant un- 
certainties. Possible models include seed BHs formed from 
Population III stars, which then gain mass through mergers 
(Schneider et al. 2002; Volonteri, Haardt, & Madau 2003) 
or through accretion (Adams, Graff, & Richstone 2001; Is- 
lam, Taylor, & Silk 2003). Alternatively, seed BHs may form 
directly from a collapsing gas cloud (Koushiappas, Bullock, 
& Dekel 2004; for a more complete discussion, see the re- 
view by van der Marel 2004). However, there are few di- 
rect observational constraints on the nature of seed BHs, the 
path by which nuclear BHs grow, or the form of the BH mass 
function. The low end of the BH mass function, in particular, 
provides critical input to check theoretical models of quasar 
formation (e.g., Haehnelt, Natarajan, & Rees 1998; Bromley, 
Somerville, & Fabian 2004), and it affects predictions of fu- 
ture gravitational wave experiments like LISA (Hughes 2002). 

A useful approach, first proposed by Soltan (1982), matches 
the local BH mass density inferred from velocity dispersion 
measurements with an accretion history inferred from obser- 
vations of luminous quasars. Yu & Tremaine (2002) conclude 
that, for standard assumptions about radiative efficiency and 
quasar lifetimes, accretion during the luminous quasar phase 
can account for the majority of the local BH mass density. 
Aside from the possibility of significant obscured accretion 
(Fabian 2004), this method is limited by its use of the SDSS 
velocity dispersion function, which at the moment is only 



well determined for early-type galaxies with er* > 200 km s 1 
(Sheth et al. 2003), corresponding to M BH > 10 8 M . 

Much of the current observational and theoretical effort on 
intermediate-mass BHs has focused on studies of star clus- 
ters (Gebhardt et al. 2002; Gerssen et al. 2002) and so-called 
ultraluminous X-ray sources (see van der Marel 2004). The 
existence of BHs in either class of objects is still far from 
proven, and, in any event, their relationship to galactic nuclei 
is unclear. With the recognition that galaxies like NGC 4395 
and POX 52 do harbor central intermediate-mass BHs, the 
most straightforward way to make progress would be to find 
many more objects like them. This study represents the first 
step in this direction. By sifting through the massive SDSS 
database, we have identified a sample of AGNs that we be- 
lieve are promising candidates for hosting BHs with masses 
below the 10 6 M threshold. 

Our sample opens up several new avenues for further study. 
While well-defined, our sample is likely to have missed 
intermediate-mass BH candidates such as NGC 4395 or POX 
52, whose spectrum is dominated by narrow emission lines, 
with the broad-line component visible only under very high 
signal-to-noise ratio. Moreover, if either AGN were embed- 
ded in a much more luminous host galaxy, the weak broad 
lines would be even harder to discern. Furthermore, we 
know that pure narrow-line (type 2) AGNs do exist, either 
because the BLR is obscured or intrinsically absent, and there 
is no reason to suppose that this situation does not persist to 
the regime of intermediate-mass BHs. A complimentary ap- 
proach, then, is to use the width of the narrow emission lines 
to select potential candidates. While the detailed geometry 
and kinematics of the narrow-line region (NLR) in AGNs are 
complex and not well understood, it has been argued that its 
velocity field is largely dominated by the gravitational poten- 
tial of the bulge of the host galaxy (Smith, Heckman, & Ming- 
worth 1990; Whittle 1992a, b). If so, then the velocities in the 
NLR are roughly virial, and on average the velocity dispersion 
of the line-emitting gas, cr g , should be approximately equal to 
the velocity dispersion of the bulge stars, ov The most de- 
tailed comparison between NLR and stellar kinematics for a 
large sample of Seyfert galaxies has been done by Nelson & 
Whittle (1996), who indeed find that in most cases <j g w <t*, 
where a g was measured using [O III] A5007. The exception 
are objects with powerful, jetlike radio structures where cr g 
can be much larger than a*, presumably because of extra non- 
gravitational motions induced by the radio outflows. Thus, in 
general one expects that <7 g , as measured in [O III] A5007 or 
other optical forbidden lines, can be taken as a rough proxy, 
or perhaps more realistically as an upper limit, for <t*. We 
are in the process of enlarging the current sample through a 
selection based on this strategy. 

Another important next step is to obtain stellar velocity 
dispersions for these objects. Such measurements will test 
whether the BH mass estimates we derive obey the M B h - <?+ 
relation. Conversely, if we assume that the M B h-c* relation is 
universally valid, the stellar velocity dispersion measurements 
will test the reliability of the virial mass estimator for AGNs. 
It would also be valuable to compare the velocity dispersions 
derived from stars with those derived from gas. Another area 
requiring attention is to obtain high-resolution, deep imaging. 
Much better data are needed to quantify the morphologies and 
detailed structural parameters of the host galaxies. Do these 
objects have bulges, and if so, what kind and where do they 
fall on galaxy scaling relations like the fundamental plane? 
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Finally, much work remains to be done, especially at non- 
optical wavelengths, to further characterize the properties of 
the AGNs themselves. 
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